Who's talking first? Consensus or lack thereof in coevolving opinion formation models 



O 

o 

(N 

> 
O 

00 



Oh! 
I 

o ■ 

o : 

o ■ 

• I— I ■ 

^ : 

Oh, 



O 



X 



Cecilia Nardini/'^ Balazs Kozma,^ and Alain Barrat-'^"'^ 

^LPT, CNRS, UMR 8627, and Univ Paris-Sud, Orsay, F-91405 (France) 
^Universitd di Padova, dipartimento di Fisica "G. Galilei" (Italy) 
^Gomplex Networks Lagrange Laboratory, IS I Foundation, Turin, Italy 
(Dated: February 2, 2008) 

We investigate different opinion formation models on adaptive network topologies. Depending 
on the dynamical process, rewiring can either (i) lead to the elimination of interactions between 
agents in different states, and accelerate the convergence to a consensus state or break the network 
in non-interacting groups or (ii) counter-intuitively, favor the existence of diverse interacting groups 
for exponentially long times. The mean-field analysis allows to elucidate the mechanisms at play. 
Strikingly, allowing the interacting agents to bear more than one opinion at the same time drastically 
changes the model's behavior and leads to fast consensus. 

PACS numbers: 89.75.-k, -87.23. Ge, 05.40.-a 



In the recent years, agent based models have been more 
and more used in the area of social sciences. Through a 
rather simple modeling approach for the individual pro- 
cesses of social influence, these models focus on the emer- 
gence of social behavior at the global population level. 
Statistical physics models and tools provide therefore a 
natural framework for such studies, and have been widely 
applied, leading to the appearance of the field called so- 
ciophysics (see [l| for a recent review on the application 
of statistical physics models to social dynamics). 

The growing field of complex networks 0, H, 0] has 
moreover allowed to obtain a better knowledge of social 
networks Bi, and in particular to show that the typical 
topology of the networks on which social agents interact 
is not regular. Many studies have therefore considered 
the evolution of models of interacting agents when agents 
are embedded on more realistic networks, and studied 
the infiuence of various complex topologies on the cor- 
responding dynamical behavior [7]. An additional fea- 
ture of networks, that may have a strong impact on the 
model's behavior, lies in their dynamical nature. They 
may indeed evolve on various timescales, and the evo- 
lution of the topology and the dynamical processes can 
drive each other with complex feedback effects. Stud- 
ies of this coevolution are more recent and still limited 
P, i, M, MMMMM,M\l^, with many open 
issues. 

In this Letter, we provide new insights into such feed- 
back effects by an investigation of Votcr-likc models 
(VM), in which agents update their opinions by imitat- 
ing their neighbors, and can also break and establish con- 
nections with other agents. More precisely, we show how 
slight modifications in the evolution rule, which have mi- 
nor consequences if the topology of interactions is kept 
fixed, can change drastically the model's behavior as soon 
as the topology can evolve on the same timescale as the 
agents' opinions. On the other hand, the simple fact of 
allowing agents to have several opinions at the same time, 
in the spirit of the Naming Game [l^ or of the AB model 
[l^, leads to more robust behavior. 

The Voter model [23| considers a population of size N 



in which each individual i has an opinion characterized 
by a binary variable Si = ±1: only two opposite opinions 
are here allowed (for example a political choice between 
two parties) psj . Starting from a random configuration 
of opinions, the dynamical evolution of the direct VM (d- 
VM) is the following: at each elementary step, an agent 
(j) is randomly selected, chooses one of its neighbors (j) 
at random and adopts its opinion, i.e. Si is set equal to 
Sj (one timestep consists of N such updates) . In the re- 
verse case (r-VM), the first agent i instead convinces its 
neighbor j (sj is set equal to Si). The distinction between 
d- and r-VM is necessary since the two interacting nodes 
do not play the same role. Moreover, the degrees of the 
first and the second chosen nodes have different distribu- 
tions, and the second is a large-degree node with larger 
probability [3]. The asymmetry in the opinion update 
between the two interacting nodes can then couple to 
the asymmetry between a randomly-chosen node and its 
randomly-chosen neighbor, leading to different dynami- 
cal properties. No important difference is expected on 
homogeneous networks but, on heterogeneous networks, 
the probability for a hub to update its state will vary 
strongly from one rule to the other. The basic imitation 
process of the VM mimics the homogenization of opin- 
ions but, since interactions are binary and random, do 
not guarantee the convergence to a uniform state. Since 
a consensus in which all individuals share the same opin- 
ion is an absorbing state of the dynamics, any finite pop- 
ulation reaches a consensus, but the time needed tc{N) 
depends on its size N and on the topology of interactions, 
and diverges as ^ oo. On static networks, tc{N) grows 
as a power-law of with an exponent depending on the 
degree distribution, and on the updating rule 21, 22, 23]. 
On homogeneous networks in particular, tc{N) oc N for 
both d- and r-VM. 

In this Letter, we consider the scenario in which agents 
can rewire their "unsatisfied" connections. More pre- 
cisely, the initial configuration is given by a random ho- 
mogeneous network of interacting agents, with average 
number of neighbors (fc) and random opinions. At each 
timestep, an agent i and one of its neighbors j are cho- 



2 




FIG. 1: (Color online) Convergence time for the r-VM as a 
function of the population size, for various rewiring probabil- 
ities. Throughout this paper, the data from simulations were 
averaged over 100 realizations of the system. Inset: same for 
the d-VM. Note the difference of scales. 



sen. With probability <&, an attempt to rewire the Hnk 
is made, if Si 4^ Sj. A new agent fc is then chosen at 
random and the Hnk is rewired to («,fc)[2^. With 
probabihty 1 — $, an opinion update takes place instead. 
The rewiring, which conserves the total number of hnks, 
is made at random: the new link is established without 
prior knowledge of the new neighbor's opinion. 

If the frequency of rewirings is small ((f> — » 0), the 
system still reaches a global consensus and the network 
remains connected. For fast rewiring rates on the other 
hand, (<i> 1) the network breaks into (typically two for 
the VM) separate connected components, each one with 
a local consensus. These two regimes are separated by a 
non-equilibrium phase transition at a critical value $c of 
the rewiring probability. Similar transitions have already 
been re port ed in co-evolving models of opinion formation 
0, [islllTj and we will not focus on this aspect here. 

A more surprising aspect of the dynamics is revealed 
by the behavior of the convergence time tc{N), which 
grows linearly with A'^ on a static network for both the d- 
and the r-VM [S^, [2^. Strikingly, the network's adaptiv- 
ity has completely opposed effects in these models (Fig. 
[1]). Consensus is strongly favored in the d-VM, for which 
tc{N) becomes cx InA^ [l^]; in contrast, for the r-VM 
tc{N) grows exponentially with the system size. The sys- 
tem therefore remains for exponentially long times in a 
state in which two groups of different opinions co-exist 
and remain connected to each other. It is noteworthy 
that the system therefore is not frozen, with agents con- 
tinuously updating their links and opinions. 

In order to understand the different behavior of the 
d- and r-VM on adaptive networks, we note that the 
state of the system is characterized by three independent 
quantities: (i) the density n+ of agents with opinion +1, 
or equivalently the magnetization m = — ri_ (n_ — 
1 — n+ is the density of agents with opinion —1); (ii) the 
number of links joining agents in the -t- opinion, Nl^-^-; 
and (iii) the number of links joining agents of opposite 



opinions, i.e. of active links Nl |_ = Nl^ (since the to- 
tal number of links is preserved, {k)/2 ~ -|-/ ). 

At the mean-field (MF) level, we can derive the evolu- 
tion equation of these quantities. Let us first consider 
the magnetization: it changes of — 2/A^ when an agent 
changes its states from -I- to — , and of +2/N in the op- 
posite case. For the d-VM, the probability of the first 
event is proportional to the density n_|- of agents in the 
-|- state, times the probability that it chooses to interact 

with a neighbor that has — opinion, i.e. /fc+ where 

is the average degree of a -I- node, and k^ = /-| /n+ 

is the average number of — neighbors of a + node. The 

probability of the second event ( > +) is obtained in 

the same way, and finally 
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In the case of the r-VM, the probabilities of the 
two processes are simply interchanged: {dm/dt)^^-^ = 
— {dm/ dt) ^yj^. On an adaptive network, it is essential 
to distinguish k+ from fc_: as shown in Fig. [21 one has 
indeed k^ > {k) > k_ if n_|_ > n_. In other words, 
the nodes of the majority opinion have more neighbors. 
This is a simple consequence of the rewiring dynam- 
ics: if TO > 0, any rewiring event — > (i,k) has 
a higher chance to randomly pick a + node as a new 
neighbor due to their larger number. Therefore, nodes 
of the larger group gain new links with larger proba- 
bility. Equation ^ then immediately shows that for 

TO > 0, {dm/dt)^YM > ^i^d ('^™/'^OrVM ^ sum- 
mary, the coevolution of opinions and topology generates 
a positive feedback for the d-VM driving the system to a 
consensus state, mutable = il; smd a negative feedback 
for the r-VM resulting in mstaUe = 0. This readily ex- 
plains the strong differences between these models. For 
the d-VM the adaptivity leads to an accelerated consen- 
sus, while it hinders the convergence for the r-VM and 
keeps the system in a dynamically evolving state with 
zero average magnetization. 

It is moreover possible to write the evolution equations 
for the various types of links. It is easy to understand 
that, according to the model's definition, the vector x — 

(rn,l-^ can evolve in 4 ways at each elementary 

update: x^x + v°, a = I,-- - ,4, with respective prob- 
abilities Let us start with the d-VM. The displace- 
ment vectors and the associated probabilities read then 
Nv^ = (2, fc„_ - /fc_+, A:_+), wi = (1 - $) n_ k_+/k_ 
Nw^ = (-2, k++-k+^, -k++), W2 = (1-$) n+fc+_/fc+ 



_/fc_; A^v^ 



Av3 = (0, - 1, 0), ii;3 

(0, —1, -t-1), Wi — '^ n\_ k^ //c+. and correspond 

to opinion changes, for which the change in magnetiza- 
tion {±2/N) is associated with changes in the densities 
of links. For example, when a — node is transformed to 

+, its links become H — and its +— links become 

++ ones (hence varies of (A;_„ — fc_| )/A^). The 

corresponding probabilities wi and W2 are obtained as 
for Eq. ([T]). v'^ and correspond to rewiring events: 



when a 



link is rewired it can be either transformed 
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FIG. 2: (Color online) Top: k+/{k) vs m for the VM. The 
symbols correspond to averages obtained from numerical sim- 
ulations with A*' = 1000, (fc) = 10, while the continuous lines 
are the numerical solution of the MF equations for the evo- 
lution of X = {m,lj^ ,'++), starting from initial conditions 

with m close to 0. Bottom: i++ and Ij, vs m. The con- 
tinuous black lines correspond to the numerical solution of 
the MF equations. The red symbols and the grey and brown 
lines correspond to single runs of the d- and r-VM, respec- 
tively {N = 500, (fc) = 10, ^ = 0.4). The inset in the middle 
shows the evolution of m for the same runs (red symbols for 
the d-VM and dashed line for the r-VM). 



to (v^) or to -l-l- (v''). For the r-VM, the displace- 
ment vectors are exactly the same as for the d-VM, but 
the transition probabilities wi and W2 are interchanged. 

and Wi remain the same since the rewiring rules are 
the same for both models. Figure [5] shows the result 
of the numerical integration of the evolution equations 
dx/dt — v"w", compared with numerical simulations 
of the models. It is clear that these equations correctly 
account for the difference between fc_|_ and fc_ and for 
the system's evolution in the phase space. Of course, the 
real systems are moreover submitted to fluctuations that 
are not taken into account in the MF description. In 
particular, looking at single runs (Fig. ^ shows clearly 
the difference between the d- and r-VM. For the d-VM, 
the density of active links decreases rapidly to and the 
system is driven to one of the consensus states. For the 
r-VM on the contrary, the system performs a random 
walk in a sort of potential well around m = with a 
non-zero density of active links, which ends only because 
of a finite-size fluctuation which leads it into one of the 
absorbing boundaries at m = ±1. 

Let us now consider that agents cannot pass directly 
from one opinion to another, but can keep both opinions 
in their "memory" , being then in an intermediate state 
that we call 0. This is the case in the Naming Game 
(NG) model, in which agents try to agree on the name 
to assign to a given object or also of the AB model 
p^ . If only two names are available (that we can call 
-|- and — for simplicity), the dynamical rules of the di- 
rect NG (d-NG) are the following: at each timestep, an 
agent i and one of its neighbors, j, are chosen at ran- 
dom to be respectively Hearer (H) and Speaker (S). S 
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FIG. 3: Convergence time for the direct (filled symbols) and 
reverse (open symbols) NG. 

proposes a name to H. If S has both names in memory, it 
chooses one at random. Let us suppose for instance that 
S proposes If H does not know the name uttered (i.e., 
it is in state — ), it absorbs this possibility by changing 
to the intermediate state, 0. If H instead has the name 
in memory (i.e., it is in state + or 0), the interaction 
is successful and both H and S agree on this particular 
name and set in state -I- after the interaction. In the re- 
verse case (r-NG), the first randomly selected agent is S 
and its neighbor is H, and the update rules remain the 
same When agents interact on a static topology, 

these dynamical rules lead to a global consensus. On 
homogeneous networks, we obtain tc{N) ^ InA^ (while 
tc{N) - N for the VM). The difference between the two 
models is due to the fact that in the VM, consensus is 
reached by a finite-size fluctuation of the average magne- 
tization while in the NG, consensus is reached due to the 
surface-tension introduced by the states, which tends 
to minimize the interface between the agents of different 
opinions and hence drive the system to a homogeneous 
consensus state [Ull^. For adaptive networks. Fig. [3] 
clearly shows that the convergence time remains logarith- 
mic for both the direct and reverse version, even if the 
r-NG is slower. The MF analysis allows to understand 
this strong difference with the VM. We can indeed write 
the evolution equation for the magnetization — n_, 
by introducing the average degree of nodes fco and the 
density of +0 and —0 links, as 




The first terms on the rhs represent the change in the 

magnetization mediated by the links. The factor 1/2 

stems from the fact that + and — nodes are not trans- 
formed instantly to their opposite counterpart but to the 
intermediate state 0. The remaining terms correspond to 
the transformation of the nodes to ± ones. For exam- 
ple, in the d-NG, the second term on the right-hand side 
is generated by the process when a node is converted to 
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FIG. 4: (dm/dt) vs m (symbols) for the r-NG. According to 
Eq. (I3J, changes in the magnetization come from r-VM-like 
interactions (dashed line) and those mediated by the 0-links, 
l+o and l-o (dash-dotted line). The upper inset gives l+o and 
;_o, the other k+/{k), k-/{k). <1> = 0.2, (k) = 10, N = lO". 

+ by first picking a node, with probability hq, then one 
of its + neighbors, with probability ko^/ko, and so on. 
Even though the first terms in Eq.-s ^ and ^ change 
sign for the d- and the r- variants of the NG just as for 
the VM, this effect is suppressed by the remaining terms 
associated with the transitions (0 — > ±) which will always 
generate a positive feedback to the change of magnetiza- 
tion. As shown in Fig. |4] indeed, /+o — '-o is of the sign 
of m, which is expected since then > This effect 
overcomes the difference between k+ and as a result, 
(dm/dt) remains of m's sign even in the r-NG, leading 
to logarithmic convergence times. Very interestingly, the 
possibility for agents to remain in an intermediate state 
before updating their opinion strongly enhances the trend 



towards consensus. 

In summary, we have shown how slight modifications 
of the interaction rules can have drastic consequences in 
the global behavior of opinion formation models in the 
case of dynamically evolving networks. In the case of the 
paradigmatic Voter model, adaptivity of the network can 
either accelerate the convergence to consensus, or on the 
contrary hinders it strongly, by maintaining the system 
in a dynamically evolving state for exponentially long 
times. A mean- field analysis allows to account for such 
differences, which are due to the coupling of the asym- 
metry between the interacting agents to the asymmetry 
in their degrees. Such coupling is known to change the 
scaling of the convergence time in heterogeneous static 
networks, which however remains a power-law of time 
^22. . l23], and in fact does not have consequences in ho- 
mogeneous networks. In strong contrast, and even if the 
adaptive network remains homogeneous, the fact that the 
majority has a slightly larger average degree suffices to 
change from a very fast convergence in logarithmic time 
for the d-VM to a dynamical state surviving for expo- 
nentially long times for the r-VM. Interestingly, if the 
agents cannot change opinion so easily, and have to go 
through an intermediate state, such as in the NG or AB 
models, convergence to consensus is enhanced also for 
adaptive networks, and irrespective of the order of inter- 
actions (d-NG vs. r-NG). The connections with nodes 
in the intermediate state determine then the dominant 
evolution of the magnetization, leading to a more robust 
behavior. 
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[26] Another possibility is that the rewired link becomes checked that this modification has no influence on the 

{j,k); the same phenomenology is then observed. overall behavior of the system described below. 
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